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Abstract 

We study Einstein-Yang-Mills equations in the presence of a grav¬ 
itating non-topological soliton field configuration consisted of a Higgs 
doublet, in Brans-Dicke and general scalar-tensor gravitational theo¬ 
ries. The results of General Relativity are reproduced in the wrd, > 
oo limit. The numerical solutions correspond to a soliton star with 
a non-abelian gauge field. We study the effects of the coupling con¬ 
stant, the frequency of the Higgs field and the Brans-Dicke field on 
the soliton parameters. 
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1 Introduction 


Einstein-Yang-Mills (EYM) equations have been investigated in several sys¬ 
tems, pp. Bartnik and McKinnon found particle-like non-abelian solutions 
of the coupled EYM theory, 121 ini- The coupling of the EYM system with a 
scalar field may lead to several theories. We mention gravitating Skyrmions, 
m, black hole solutions in dilaton, 0^1, and massive dilaton and axion grav¬ 
ity, 13, and other field configurations in the EYM-Higgs theory with a Higgs 
doublet, BM or with a Higgs triplet, Iiniiiiiii2i- 

Non-topological solitons and soliton stars have a long presence in modern 
physics, uni- In the present work we need the theory of q-balls and q-stars. 
Q-balls are non-topological solitons in Lagrangians with a global U{1) sym- 
metry, jT7j, or an SU{3) or SO{3) symmetry, jTH]. Their relativistic general¬ 
izations may consist of one or two scalar fields, I2DI. in a Lagrangian with a 
global U{1) symmetry, or of a non-abelian scalar field in the adjoint represen¬ 
tation of SU{3), |2II, or of a scalar and a fermion field, [22] in asymptotically 
flat or anti de Sitter spacetime, j2S|- Q-solitons with local symmetries have 
also been investigated. There are charged q-balls, m, and charged q-stars, 

|21j. 

Interesting alternative gravitational theories are the scalar-tensor gravi¬ 
tational theories, which appeared in the original paper of Brans and Dicke 
[12], where the Newtonian constant G was replaced by a scalar held 0 bd, 
and the total action contained kinetic terms for the new held times an cubd 
quantity. (Ubd was regarded as a constant in the original paper. The theory 
generalized in a series of papers, 131121, mainly in the direction of replac¬ 
ing the constant uj-qd with a function of the Brans-Dicke (BD) scalar held. 
Within the BD gravitational framework, Gunderson and Jensen investigated 
the coupling of a scalar held with quartic self-interactions with the metric 
and the BD scalar held, 0BD, |22|- The properties of boson stars within this 
framework have been extensively studied in a series of papers pHl 1771 |2H| • 
Their results generalized in scalar-tensor gravitational theories, 1221 EOj- The 
case of charged boson-stars in a scalar-tensor gravitational theory has been 
analyzed in HQ. 

The purpose of the present work is to hnd numerical solutions, resembling 
q-stars, of the EYM equations in the presence of a Higgs doublet in the fun¬ 
damental representation of SU{2) in BD theory or in a general scalar-tensor 
gravitational theory and to compare the above results with the solutions ob¬ 
tained in General Relativity. In the absence of the gauge held, the equations 
of motion give rise to a gravitating non-topological soliton, when using a cer- 
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tain potential for which oje = \/t//|*hp ■ < m, where m is the mass of the 
free particles, imposing an harmonic time-dependence on the scalar held and 
equalizing the frequency to Our gravitating soliton is non-topological in 
the sense that <h, t/ —>• 0 for p —>• cx) according to ng. It is a q-type non- 
topological soliton in the sense that in the absence of both gravitational and 
gauge helds one can hnd by simple calculations that this spherically symmet¬ 
ric Higgs held rotates within its symmetry space with a frequency oje equal to 
the minimum of the quantity, as in q-balls. The diherence between 

this soliton and the usual non-abelian q-balls is that the symmetry space in 
the case of non-abelian q-balls is the entire SU{3) space but in our case is an 
abelian U{1) subgroup of the SU{2) group, though both held conhgurations 
are non-abelian. Also, we investigate the fundamental and not the adjoint 
representation of SU{2). In any case we hnd an analytical solution for the 
scalar held within the soliton, using the approximation known by the study 
of q-stars. 

2 EYM equations in Brans-Dicke gravitational 
theory 

We consider a Brans-Dicke scalar held 0 bd, a matter Higgs scalar doublet, 

$, in the fundamental representation of SU (2) and a gauge held A, coupled 
to the metric The total action is: 

S — J ^ ^v~9 ( 0BD-R — ^- 1 

T d X-\/ p/Imatter ) (1) 

with (Ubd a constant in BD gravity and a certain function of the 0 bd held in 
generalized scalar-tensor gravitational theories and: 

-^matter = ^TrF^.F'^" + (D^<I>)t -U . (2) 

In the above Lagrangian we dehne: 

- dyA^, - i[A^, A^] . 

One may use the = d^Ai, — d^A^ — «g[A^, A^] form and reproduce eq. 

El rescaling: —>• A^/g with g the gauge coupling, or held strength. The 
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one-form gauge field A is: A = A^dx^ = TaA’^^dx^, with = |ra and Ta 
the Pauli matrices. The factor K appearing in the action is dehned by the 
relation Tr(TaTb) = K6ab, reading K = 1/2. 

In the presence of the Brans-Dicke scalar, the Einstein equations take the 
form: 


G 






g^J.u4>BD■x 




CUrd 



di^4>BDdu4>BD — 


(4) 


and the equation of motion for the BD field is: 


2cc;bd 



;A 


;A d^4>BI)d\(l)BD . n 

-OJbD -To- R 

0BD 


0 . 


(5) 


is the Einstein tensor, R is the scalar curvature and T^,y the energy mo¬ 
mentum tensor for the matter helds (gauge and Higgs) given by the equation: 

T,, = ^Tr + 

- g,,U . (6) 


Tracing Einstein equations and substituting the result in the Lagrange equa¬ 
tion for the BD held we take: 


;A 


Stt 


2a;BD + 3 


-T . 


(7) 


We will choose a general, spherically symmetric held conhguration, dehn- 
ing: = (sin'd cos 99 , slu'd sin 99 , cosd) and Tp = rATa, = i9,?Tp and 

T<p = (1/sin'd)(9(pTp. The gauge held and the Higgs doublet are: 


A = aTp -(- i{l — Re(x;)[Tp, dTp] -|- Imu;Tp = 

aTp -|- [Ima;T.ij -|- (Recu — l)T(p]d'd -1- [ImcuTip 4- (1 — Reci;)T,j] sin-dd^?, (8) 


$ = o-exp(*^Tp)|6) (9) 

with a = a{p,t), ^ \b) a constant unit vector of the internal SU{2) 

space of the scalar (Higgs) held and a = a^dt + Qpdp. In order to form 
a held conhguration corresponding to a charged q-type soliton we choose: 
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flo = oo(p), ct(p, t) = a{p) and ^ = UEt, when we choose = 0 for simplicity. 
The ansatz a{p,t) = a{r) and ^ = UEt is the obvious generalization to the 
(j){p, t) = (j(p)e*^* ansatz, known from q-solitons. So, the role of the eigen- 
frequency, n is now played hy ue- The choice Op = 0 implies that our field 
conhguration is not very general, but our purpose is not to hnd the more 
general solution, but a proper one, with the above features and resulting to 
stable solitons. The BD scalar held is supposed to be static and to posses 
a spherical symmetry, as the matter held conhguration. With the above 
assumptions, we can write a static, spherically symmetric metric: 

ds"^ = + p^ sin^ . (10) 


The matter action takes now the form: 




matter — 


p^ sin d 

~7W L 2g: 


a'^AB + 2 


\^\ % ^ (kr -1)^ 


p^ 


p^ 




1 ^ 

■-{ue - ao^a'^B + ^^[(Reci; - cos{uEt)Y + (Imu; - sin(a;£;t))^] - U 
4 2p^ 


( 11 ) 


In order the action to be time-independent we may choose Recu = cos{uEt) 
and Imu; = sin(a;£;t), but this choice is not a solution to the equation of 
motion for a;, or ct; = 0 which is a solution of the equation of motion, so, our 
solution is embedded abelian. 

If m is the mass of the free particles, we make the following rescalings: 
p = 2mp, = a, = — , a=^, 

p = , g = g€“‘ , 4 bd = + , (12) 

2u;bd + 3 

with: 

e = V^TiGrri? . (13) 

We dehne: 



Gravity becomes important when R ~ StiGM, with < 0 >= 1/G. With 
our rescalings and because the energy density within the soliton is ~ m^, we 
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find that r ~ 1 and if a varies very slowly within the soliton, from a (t( 0) 
value at f = 0 to a zero value at the outer edge of the soliton surface, then 
V ~ e^m^. For m ~ GeV the 0(e) quantities are negligible. We choose a 
simple rescaled potential, admitting q-ball type solutions in the absence of 
gravity and gauge helds, namely: 

U = , (15) 

V 3m^ / 


which with our rescalings and after some algebra takes the form: 



(16) 


where we set m = 1. From now on we drop the tildes and the 0(e) quantities. 
From the equation of motion for the Higgs held we hnd: 


a^ = l+eoB^/\ t/ =1(1 +0353/2) ^ 

T = 2hF - 40 , 


W = 0^B(1 + , 


(17) 


with: 

^0 = — Oo • (18) 

The equation of motion for the Higgs held within the surface gives a 
boundary condition for the gauge held 6*o, which reduces to an eigenvalue 
equation for the frequency in the case of global SU{2) symmetry (i.e.: when 
ao = 0). The surface width is of 0(m“^). The Higgs held a varies rapidly 
from a cxo value at the inner edge of the surface, to a zero value at the outer 
one. Dropping form the Lagrange equation the 0(e) terms and integrating 
the resulting equation, we hnd that within the surface: 


V + W-U = 0 . 


(19) 


In order to match the interior with the surface solution we set cr' = 0 at the 
inner edge of the surface. Then, using eqs. [T7| and UHl we hnd: 


^Osur — 


. 1/2 0 - 1/2 
-OLsur -^sur 


( 20 ) 


where ^osur is the value of 6o within the thin surface. In the absence of gauge 

1 /9 

helds we take: uje = Asur/2 which in the absence of gravity gives uje = 1/2, 
which is the correct eigenvalue equation for the q-soliton frequency. 
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With the above definitions, the independent Einstein equations take the 
following form: 


A-l IdA 2o;bd + 3 / 2W - AU\ 

r dr (2a;BD + 4)$BD v 2a;BD + 3 / 

/d<hBD\ _ A dB d^BD 

2<h|j3 \ dr ) 2 ^bdB dr dr 


( 21 ) 


A — 1 A1 dB 2ci;bd “1“ 3 /rir tt 2W — AU\ 

-=- { W — U - 

r 2 B r dr (2a;BD + 4 )<I'bd V 2a;BD + 3 / 

CUBD^ / d$BD \ ^ A f d^4>BD 1 dy4d$BD \ ^22) 

\ dr J <I>bd V dr^ 2A dr dr J ' 


the Euler-Lagrange equation for the BD scalar is: 

^ r d^$BD /2 dd>BD l ^ 2W-4U 

dr‘^ \r 2A dr 2B dr) dr 2 uibb) + 4 ’ 


and the equation of motion for the new gauge field 6q\ 



[2 ^ A\ B’\„, gX(l + flo-S"''") 

\r'^ 2A^ 2 b) ” 2^ 


0 . 


(24) 


with boundary conditions: 


24(0) = !, 0'(O) = O, e{R) = e,^,, 4>bd = 1, (25) 


and for r —>■ oo\ 


A{r) = l/B{r) = 1 , ^bd = 1 , Oij) = uje ■ (26) 

The first condition of eq. [^reflects our freedom to redefine 24(r), the second 
and forth result from the spherical symmetry of the configuration and the 
third is the eigenvalue equation for the new gauge field Oq, with R the soliton 
radius. The first condition of eq. |2S1 is a straightforward consequence of 
the Einstein equations for localized matter configurations, the second is the 
boundary condition for the BD field with the proper rescalings and the third 
denotes the absence of gauge fields at infinity. We numerically solve the 
coupled system of eqs. I^TII^ 
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Figure 1: The value of the Higgs held a at the center of the soliton as a 
function of the coupling constant g^. The numbers within the hgures fTlIKl 
denote the eigenvalue ^osur, which for g^ = 0 reduces to a;^. Dashed lines 
correspond to (Ubd = 5 and solid lines to cubd = 500. The results of General 
Relativity almost coincide with the BD theory for cubd = 500. 



Figure 2: The radius of the soliton as a function of g^. 








Figure 5: The value of the gauge held, oq, at the center of the soliton as a 
function of 

The Noether currents corresponding to the generators of the SU (2) alge¬ 
bra are given by the relation: 


_ dL dL 

JOa I d{dQ^) d{dQ^*) 


iTa 0 
0 —iTr, 


$ 

<|)* 


(27) 


We can hnd that: 


joi = ^a'^Oosm^cosip , 
302 = ^n^6'osin?9sin(^ , 
Jos = ^ 0 -^610 cosh* , 


(28) 


and 


Jo = \lill + Jo 2 + Jos - ■ 


The particle number is: 


(29) 


N = 2tt f a'^Ooy ^r'^dr . 


(30) 


The total energy of the held conhguration results from the relation: 


,, , ^ 2GM , Gg^N^ 

A{p) = 1-+ ——— , p^oo , 


P 


Anp^ 


( 31 ) 
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which gives with our rescalings: 


/ g^iV^ \ 

M = 47rr M - A{r) + 1 , r ^ oo . (32) 

For cubd oo the Einstein equations take the simple form: 


A - 1 

A' 

n/2 

-U-W- ^AB , 

2g2 

(33) 


r 

A-1 

AB' 

ni2 

= W-U-^AB , 

2g2 

(34) 

^2 

B r 


when the other relations remain the same and the results of General Rela¬ 
tivity are reproduced. 


3 EYM equations in general scalar-tensor the¬ 
ory 


In the original BD gravitational theory cubd is a constant. In a more general 
theory it may be regarded as a function, usually of the BD held. We will use 
one of the forms investigated in a cosmological framework, isainn], namely: 

2a;BD + 3 = (Uq^bd ) (35) 


with ujq and n constants. The Lagrange equation for the BD held is: 


;A 


UjQCf) 


n 

BD 


SttT 


du 


BD 


d(f>BD 


>P 


We rescale: 


Cl>o — 


( 


2ujbd + 3 




(36) 


(37) 


V 2a;BD + 4 ^ 

and the other quantities as in[T21and drop the tildes and the 0(e) quanti¬ 
ties. We will use n = 1. For other values of n the behavior of the soliton 
parameters is very similar. The Einstein equations take the following form: 


R- 1 


1 dA 
r dr 


UJq 


<a^o4*bd + 1 
2W -AU 


-W -U 


CUo^bD 
^4>bd lA;o‘hBD + 1 


X 




BD 


(^o*hBD — 3 R’I’bd AB'A? 


■ BD 


2<1>Id 


2^bdB 


(38) 
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Figure 6: The value of the held a at the center of the soliton as a function 
of uq. The numbers within the hgures IHlfTUl denote the g^ value. We use 
6*osur = 0.4. In hgures Emni solid lines are the numerical results produced in 
the framework of scalar-tensor gravitational theory, when dashed lines are 
the results from General Relativity 



Figure 7: The radius of the soliton as a function of uoq 
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Figure 8: The asymptotically anti de Sitter mass, M of the soliton as a 
function of ooq. 



Figure 9: The particle number of the soliton as a function of ujq. 
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Figure 10: The value of the gauge held at the center of the soliton as a 
function of ujq. 


A - 1 AldB 


Uq 


B r dr cuo^bd + 1 


2W-AU 


W-U 


CUo^bD 
y4$BD O^o'^BD + 1 


X 


$ 


BD 


o;o4>bd - 3 24<FgD ^ ^ 


2^Id 


4> 


BD 


24> 


BD 


The equation of motion for the BD held is: 

A 


Kb + 


2 ' 
2A 2b) 


1 


O.'O^BD + 1 


2W -4U- 


^“^BD i^O^BD + 1 


4> 


BD 


(39) 


(40) 


The equations of motion for the gauge and Higgs held remain unchanged. 
We solve numerically the coupled system of eqs. IHTIHHI and |211 when the 
relations EDI IS21 for the soliton mass and particle number hold true. 


4 Conclusions 

We studied EYM equations in the presence of a Higgs doublet in Brans-Dicke 
and a simple scalar-tensor gravitational theory and compared our results with 
the solutions of General Relativity. The Higgs doublet is characterized by a 
potential admitting q-star type and and q-ball type solutions in the absence 
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of gauge fields and gravity and gauge fields respectively. It is a matter of 
simple algebra to verify this claim. So, the EYM-Higgs equations reduce to a 
system of equations corresponding to (charged) soliton stars. These objects 
are stars and not black holes, having no horizon or other anomalies. 

There are two crucial parameters, resulting from the soliton star itself, 
the held strength, g and the eigenvalue 6'osur which in the absence of gauge 
helds reduces to the usual soliton eigen-frequency. The above eigenvalue is 
straightforward connected to gravity strength on the surface, through eq. 1201 
So, a soliton with small 6*osur shows a stronger gravitational force on its sur¬ 
face, which corresponds to a more massive, or denser soliton, and this can be 
verihed by our hgures. Also, larger value for the held strength, increases the 
energy and radius of the held conhguration due to the electrostatic repulsion 
between the diherent parts of the soliton, when the value of the scalar held 
within the soliton decreases for the same reason. These results hold true in 
General Relativity as well as in scalar-tensor theories. 

In hgures ITIIHI the results of General Relativity are not depicted because 
they almost coincide with the cubd = 500 (solid lines) case, cubd — 500 is 
the lower experimental limit. The results of general relativity are exactly 
reproduced when cuBD —>• oo. In hgures IHlITni we study the behavior of the 
soliton parameters for 5 < ujq < 1000. The results of General Relativity are 
practically reproduced for uiq = 1000 as one can see from the dashed lines in 
hgures IHHTHl which correspond to the results of General Relativity. 
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